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Abstract 



Airless planetary bodies are covered by a dusty layer called regolith. The grain size of the regolith determines the temperature and 
the mechanical strength of the surface layers. Thus, knowledge of the grain size of planetary regolith helps to prepare future landing 
and/or sample-return missions. In this work, we present a method to determine the grain size of planetary regolith by using remote 
measurements of the thermal inertia. We found that small bodies in the Solar System (diameter less than ~ 100km) are covered 



by relatively coarse regolith grains with typical particle sizes in the millimeter to centimeter regime, whereas large objects possess 
very fine regohth with grain sizes between 10 and 100 jum. 
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1. Introduction 



Planetary surfaces are exposed to a continuous flux of im- 
pactors of various sizes, which have, due to the hyper-velocity 



nature of the impacts, ground down the initially rocky mate- 
rial to ever finer particle sizes ([Chapman 



1976 Housen and 



P- |Wilkeningl [1982) 1. Thus, airless bodies in the Solar System 



Q (e.g.. Mercury, the Moon, planetary moons, and asteroids) 
Yi are covered by planetary regolith, a granular material consist- 
ed ing of distinct solid grains (Chapman 2004 1. Each 
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hyper- 
velocity impact causes the ejection of material from the forming 
crater, which can partly be re-captured by the planetary body 
if its gravitational escape speed exceeds the ejection velocity 
( [Housen and Wilkening"| [1982[ [Cintala et aL| [1978| l . Fragments 
with higher velocities are lost into interplanetary space. 
Laboratory experiments ([Fuj 



Nakamura and Fujiwara 1991 



iwara and Tsukamoto! '19801 
Nakamura, 1993; Nakamural 



et al.[|1994[ l and studies of the crater structures on the Moon 



(Vickery 1986 1987|l have shown that hyper- velocity impacts 



accelerate small fragments to higher velocities than large ones 
so that we expect that larger planetary bodies possess a finer av- 
^ erage regolith grain size than small bodies ( jCintala et al.[|1978] l. 
' However, the mean particle size of planetary regolith has, ex- 
cept for the Moon ( jPuke et aL 



X 



1970 McKay et al. 2009 1, not 



been measured directly. 

In this work, we present a new method to determine the grain 
size of planetary regolith from remote or in-situ thermal inertia 
measurements (see Sect. [2] and Fig. [T]i. Using literature data 
of the heat conductivity of lunar regolith ( [Cremers et al.[ [1971[ 
Cremers and Birkebak[ [1971[ ), sampled by the Apollo 1 1 and 



12 astronauts, we fit our published heat-conductivity model of 
granular material ( Gundlach and Blum 2012 see Sect. [3]) to the 
lunar regohth and show that the resulting characteristic grain 



size agrees with the regolith size distribution measured for the 
Moon (McKay et al.[ [2009[ see Sect. |4|. We confirm the ex- 
pected anti-correlation between the regolith grain size and the 
gravitational acceleration of the planetary body for a large num- 
ber of asteroids, the Moon, the Martian moons and Mercury 
with diameters between 0.3 km and 4,880 km (see Sects. [5]|8]l. 
Our results can help to prepare future landing and sample-return 
missions to primitive bodies of the Solar System (see Sect. |9|. 

2. Strategy 

To determine the particle size of planetary regolith (see Fig. 
[T]for a flow chart of our approach), we used literature data of 
thermal inertiae measured for different airless bodies in the So- 
lar System (see Sect. |5]). The thermal inertia. 



(1) 
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describes the resistance of the near surface material of a Solar 
System body to follow diurnal changes in the irradiation and 
depends on the heat conductivity A and the volumetric heat ca- 
pacity of the bulk regolith, respectively. The volumetric heat 
capacity can be expressed by C = 0p c, with the packing frac- 
tion 0, the mass density p, and the specific heat capacity c of 
the regolith particles. Thus, the thermal inertia depends on the 
material properties of the regolith and the degree of compaction 
of the regolith particle layers. We also expect that the heat con- 
ductivity depends on the size of the regolith particles as sug- 
gested by [Gundlach and Blum ( 2012[ l. With the knowledge of 
the thermal inertia, the surface-material properties and assump- 
tions about the packing density of the regolith particles, one can 
thus derive the particle size of the surface regolith. 

For the estimation of the mass density and the specific heat 
capacity of the regolith material, we divided the analyzed bod- 
ies into stony (S), carbonaceous (C), and metallic (M) objects 
(see Sect. |5]l. For each of the three classes, the mass density and 
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Figure 1 : Strategy for the grain size determination of planetary regolith using 
thermal-inertia measurements. 



heat capacity of the material were approximated by the labora- 
tory measurements of these properties of representative mete- 
orites ( ^Opeil et al. 2010| , i.e. Cronstad and Lumpkin (S), Cold 



Bokkeveld and NWA 5515 (C) and Campo del Cielo (M) (see 
Sect. [3]). In order to take different degrees of compaction of the 
regolith into account, we treated the volume filling factor (or 
packing fraction) of the material as a free parameter and varied 
it between = 0.1 and - 0.6 in intervals of A(;* - 0.1. For 
each of the six packing densities, we thus arrive at a heat con- 
ductivity value. In general, the heat capacity and the thermal 
conductivity are temperature dependent. Thus, we also derived 
the surface temperatures of the celestial bodies in our sample 
at the time of observation of their thermal inertiae. The mean 
grain size of the planetary regolith can then be determined from 
a comparison with a modeled heat conductivity of granular ma- 
terials in vacuum ( |Gundlach and Bluml|2012| see Sect. [6]l. 

3. The heat conductivity model 

The heat conductivity model for regoUth in vacuum used in 
this work was introduced by'Gundlach and Blum|(2012l and 



is an extension of the model by Chan and Tien (197 3)1. Re- 
golith is a packing of individual solid particles, in which heat 
can be transported through the solid network of grains and due 
to radiation inside the pores of the material. Due to the rel- 
atively small contact areas between neighboring regolith parti- 
cles, the heat conductivity through the particle network is gener- 
ally much lower than the heat conductivity of the solid material 
the particles consist of. Radiative heat conduction is favored by 
large void spaces between the regolith particles, i.e. by a large 
mean free path of the photons. 

The heat conductivity of a granular packing of equal-sized 
spheres with radii r temperature T and packing fraction is 
given by 

A(r, T, 0) = A,om(T) H(r, T,cp) + % tr eT^ A(r, 0) (2) 



and consists of a conductive (first term on the rhs of Eq. [2| 
and a radiative term (second term on the rhs of Eq. |2|. It re- 
lates the granular heat conductivity to that of the solid material 
of the regolith particles, /Isoiid through the dimensionless Hertz 
factor H{r, T, (p) (describing the reduced heat flux through the 
contacts between the regolith particles) and the mean free path 
Mr, (p) of the photons within the pore space of the granular ma- 
terial. Here, cr and e denote the Stefan-Boltzmann constant and 
the material emissivity of the regolith grains, respectively. The 
mean free path of the photons can be expressed by 



1 - (j) 
A(r,(f>) = ei — - — r. 



(3) 



wifliei = 1.34 ± 0.01 ( |Dumen| [TWT] [Gundlach and Blum| 
2012| l. Thus, the radiative heat conductivity is proportional to 



the regolith grain size. 

Due to the expected smallness of the regolith particles and the 
low gravity environment for asteroids, the inter-particle forces 
in the uppermost regolith layers are dominated by mutual van 
der Waals attraction between the regolith grains and not by the 
particle weight. Thus, the Hertzian dilution factor for such a 



granular packing can be expressed by (Gundlach and Blum 
[2012) 1 



Hir,T,<p) = 



971 1 - y{T) 

~4 ~E ~ 



1/3 



ifiexp[f2<p])-X- (4) 



The first of the three dimensionless factors on the rhs of Eq. |4] 
describes the heat-flux reduction due to the neck between the 
particles caused by the van der Waals force. Here, yu, E, and 
yiT) are Poisson's ratio. Young's modulus, and the specific sur- 
face energy of the grain material, respectively. The derivation 
of this first dimensionless factor is based on the solution of the 



heat transfer equation for solid spheres in contact. [Chan and 



Tien ( 1973| l extended this to a network of monodisperse spheres 



in contact. The resulting dilution factor depends on the external 
applied force acting on the particles. This force determines the 
Hertzian contact area between the particles and, hence the effi- 
ciency of the heat conduction through the network of particles. 
For large particles on Earth (r ~ 50 fim), the heat conductiv- 
ity can be derived using the weight of the particles. However, 
for small particles (r ~ 50yL(m) on Earth, the adhesion force 
(e.g., van der Waals bonding) between the particles is much 
larger than their weight. For asteroids, for which the gravita- 
tional force is comparatively small, adhesion between the par- 
ticles is always the dominant force in the first few regolith par- 
ticle layers, determining the contact area between the regolith 
particles. In order to implement the adhesive force between 



the particles into the model, we used the JKR theory (Johnson 



et al.|[ 197 11 for the calculation of the force between particles in 



contact. Using the model by |Chan and Tien ( 1973) 1 and replac- 
ing the external force by the internal JKR force, yields the first 
dimensionless factor on the rhs of Eq. [4| 

The second factor on the rhs of Eq. [4[ contains the structural 
information about the particle chains inside the regolith, which 
can transport heat, with /i = (5.18 + 3.45 ) x 10"^ and /2 = 



2 



Table 1 : Summary of the material properties used in the model. 



Physical Propeity 



Value 



Unit 



Ref. 



Material 



E 

y(T) 

'Isolid.S 

'lsolid,c(7') 

'lsolid,M(7') 

Ps 

PC 
PM 

es(r = 200 K) 
cciT = 200 K) 
eM(r = 200 K) 
es.c 



0.25 

(7.8 ± 1.9) X 10'» 
6.67 X 10-^T[K] 
2.18 

1.19 + 2.1 X 10-^ r[K] 
12.63 + 5.1 X 10"- T[K] 
3700 
3110 
7800 

500 

560 

375 

1.00 

0.66 



[Pa] 

[Wm-' K-'] 
[Wm"' K"'] 
[kgm-3] 
[kgm--^] 
[kgm--'] 
[Jkg-'K-i] 
[Jkg-'K-'] 
[Jkg-'K-i] 



Schultz 



Schultz 



1995 



1995 



Gundlach and Blum 2012 



Opeil et al. 


2010 


Opeil et al. 


2010 


Opeil et al. 


2010 


Opeil et al. 


2010 


Opeil et al. 


2010 


Opeil et al. 


2010 


Opeil et al. 


2010 


Opeil et al. 


2010 


Opeil et al. 


2010 



Basalt 
Basalt 
Si02 

Cronstad and Lumpkin (Meteorites) 
Cold Bokkeveld andNWA 5515 (Meteorites) 
Campo del Cielo (Meteorite) 
Cronstad and Lumpkin (Meteorites) 
Cold Bokkeveld and NWA 5515 (Meteorites) 
Campo del Cielo (Meteorite) 
Cronstad and Lumpkin (Meteorites) 
Cold Bokkeveld and NWA 5515 (Meteorites) 
Campo del Cielo (Meteorite) 



/i: Poisson's ratio. 

E: Young's modulus. 

y(T): specific surface energy. 

/lsoiid(7'): heat conductivity of the solid material. 

p: density of the solid material. 

c: heat capacity of the solid material. 

e: emissivity of the material. 



Physical Property 



Table 2: Material properties used for the calibration of the heat conductivity model with lunar regolith. 



Value 



Unit 



Ref. 



Material 



r 


20 




f 


0.25 




E 


(7.8 ± 1.9) X 10'" 


[Pa] 


y(T) 


6.67 X K)-' r[K] 


[Jm-2] 




1.78 


[Wm-' K- 


Pregolith 


1300 


[kgm-3] 


PsoHd 


3040 


[kgm-3] 


<P 


0.43 





McKay etal. 2009 



Schultz 



Schultz 



1995 



1995 



Gundlach and Blum 



2012 



Cremers et al. 


19 


71 


Cremers a 


-id Birkebak 


1971 




Warren||2001 





Lunar Regolith 
Basalt 
Basalt 
Si02 

average between SIOt and meteorites 
Lunar Regolith 
Lunai" Meteorites 



1 



r. mean radius of the grains. 

fi: Poisson's ratio. 

E: Young's modulus. 

y(Ty. specific surface energy. 

^soiidiT): heat conductivity of the solid material. 

Pregolith^ density of the regolith. 

Psoiid^ density of the solid material. 

0: volume filling factor of the regolith. 

6: emissivity of the material. 



3 



0.94 being empirical constants (Gundlach and Blum 



5.26 ± 
|20T2] l. 

Finally, the factor x describes the reduction of the heat con- 
ductivity when the model assumption of monodisperse spheri- 
cal particles is relaxed towards irregular polydisperse grains of 
a real regolith. This factor is determined by the calibration of 
the heat conductivity model using lunar regolith (see Sect. |4]). 

Using Eqs. |2]-|4] the heat conductivity of regolith in vacuum 
can be expressed as 



A(r, T, 



Isolid 



9n 1 - y(T) 
~4 ~E 7~ 



1/3 



x(/iexp[/2 0])-;r + ScreT^ei 



(5) 



Eq. [5] shows the dependence of the heat conductivity on the 
grain size of the regolith, which is the only free parameter of the 
model if the material properties, the regolith temperature, and 
the volume filling factor of the regolith are known. For small 
regolith particle sizes, the first term on the right hand side of Eq. 
[sjdominates so that /I oc r '^^. For large regolith particles, high 
temperatures, and/or very high regolith porosities, the second 
term dominates, so that we expect A cc r. 

The material properties used to evaluate Eq. |5] are summa- 
rized in Tab. [T] For Poisson's ratio and Young's modulus, we 
used the values measured for basaltic rock. The bulk heat con- 
ductivity was estimated from measurements of the heat con- 
ductivity of meteorites ( Opeil et aT] |2010| ). Therefore, the ana- 
lyzed objects were divided into three different spectral classes, 
S (stony objects), C (carbonaceous objects) and M (metallic 
objects). For each class, representative meteorites were cho- 
sen, i.e. Cronstad and Lumpkin for class S, Cold Bokkeveld 
and NWA 5515 for class C, and Campo del Cielo for class M, 
respectively. Since the meteorite samples are porous, we cal- 
culated the bulk heat conductivity of the non-porous material 
using Maxwell's formula ( [Maxwell and Thompson! |1904| |, 

2 + 1// 



■*^s olid, non- porous 



= A. 



solid.porous 



2(1 - ^) 



(6) 



Here, iff - 1 - is the porosity of the respective meteorite 
(i/r = 0.167 for Cronstad (S), ^ = 0.191 for Lumpkin (S), 
i// = 0.373 for Cold Bokkeveld (C), if/ = 0. 251 for NWA 5515 
(C) and (A = 0.012 for Campo del Cielo (M), |Opeil et aL]|2010l l. 
Maxwell's formula was derived for a solid body containing a di- 
lute "suspension" of voids, which is a more appropriate descrip- 
tion for the structure of a meteorite samples than the regolith 
model. Eq. |6] together with the measured heat conductivities 
of the meteorites from 'Opeil et al. (2010 1, yield the bulk heat 
conductivities of the non-porous materials used in this work. 



A 



solid, S 



2.18Wm"'K"^ 



AsoUd,c(T) = (1.19 + 2.1 X 10-^' r[K]) Wm-' K"' 



A 



solid, M 



1 t'-i 



(T) = (12.63 + 5.1 X 10"^r[K])Wm"'K 



(7) 



Since the temperature dependence of the specific surface en- 
ergy of basaltic rocks and meteorites is not known, we used the 



corresponding value estimated for Si02 (Gundlach and Blum 
2012[ l. Furthermore, an emissivity of e = 1 was assumed for 
the S and C class asteroids. The calculations for the M class 
asteroids were performed with an emissivity of e = 0.66. This 
value was estimated using Eq. |8]together with the measured sur- 
face temperature and the albedo of the M class Asteroid Steins 
(geometric albedo: 0.22; Bond albedo = 0.10; see Sect. |5]for 
details). 



4. Calibration of the heat conductivity model using lunar 
regolith 

In order to take the irregular shape of the regolith particles 
into account, we determined the parameter in Eq. [4]and[5]us- 
ing ground-based measurements of the heat conductivity of lu- 
nar regolith and, thus, re-calibrated the heat conductivity model 



( Gundlach and Blum 2012j l to natural regolith. 

For the caUbration of the heat conductivity model, we used 
the measured temperature dependence of the heat conductivity 
of lunar regolith returned with the Apollo- 1 1 and the Apollo- 12 
missions (Cre mers et al.| |1971[ [Cremers and Birkebak| |1971| l 
and compared it to the result of our heat conductivity model 
(see Sect. [3]), using the known or estimated material proper- 
ties of the lunar regolith samples (see Tab. |2]). For complete- 
ness, we should mention that the heat conductivity of lunar re- 
golith was also estimated from radio emission measurements of 
the lunar surface ( |Krotikov and Troitskii[[T963l l. However, the 
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Figure 2: Calibration of the lieat conductivity model witfi lunar regolith. The 
measured heat conductivity of lunar regolith is plotted as crosses (ApoUo-ll; 
[Cremers et al.| [T971 1 and pluses (Apollo-12; [Cremers and Birkebak] |1971^ . 
Our heat conductivity model of granular material (Eqs. |2|-|5J was fit to the 
data in order to estimate the model parameter^, which takes the irregular shape 
of the particles into account. The solid and the dashed curves show the model 
results for x = 0.41 without and with radiative heat transport, respectively. 
For comparison, the model results foTX = ".31 and^ = 0.51 are also shown 
as dotted curves. For these calculations, a volume filling factor of </> = 0.43 
jCremers et al. 1971[[Cre mers and Birke bak[[197H and a mean particle radius 



20 ^m (McKay etal. , , 2009^ 1 were used. Additionally, the dashed-dotted 



of r 

curve shows the radiative heat conductivity without the network heat conduc- 
tivity for comparison. 
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heat conductivity derived from radio emission measurements 
is approximately one magnitude higher compared to the val- 
ues measured in the laboratory. Because of the better accuracy 
of the laboratory experiments (the error of the radio emission 
measurements was estimated to be 50 % of the measured value, 
IKrotikov and Troitskii 1963 ), we decided to use the results ob- 



tained byjCremers et al.,(,1971j and the results obtained by Cre- 
Imers and Birkebak ( 1971 1. 



Fig. |2] shows the temperature dependence of the measured 
heat conductivity of lunar regolith returned with the Apollo- 1 1 



(crosses; [Cremers et al. |1971 1 and Apollo-12 missions (pluses; 



Cremers and Birkebak 197 1 1 together with the result of the heat 



conductivity model described in Eqs. 2]|5 without (solid curve) 
and with radiative heat transport (dashed and dotted curves). 
The model parameters used here are summarized in Tab. [2] 

The volume filling factor (f) - 0.43 was estimated us- 
ing the meas ured density of the lunar regohth, Pregoiith = 



1300 kg m ^ (Cremers et al. 1971 Cremers and Birkebak 



[1971 ), and the assumption that the density of the solid re- 
gohth material is represented b y the density o f the lunar me- 



teorites, psoiid - 3040 kg m""* (Warren 2001 



For the heat 

conductivity of the solid mater ial, we used the two limiting val- 
ues isoiid = 2.18Wm'' K"' (|Opeil et al.| |2010| and A^oViA = 
1.37 Wm ' K ' for SiOa (|Ratcliffe| |1963]l. We obtain for the 



only unknown parameter in the model, which describes the 
reduction of the heat conductivity when spherical grains are re- 
placed by irregular grains, values of x - 0-31 ± 0.01 x - 
0.50+0.02, respectively. Thus, we estimate a mean thermal con- 
ductivity of the regolith material of /Isoiid = 1 .78 W m"' K"' and 
a mean reduction factor of = 0.41 ± 0.02 (stat) + 0.10 (syst). 
Fig. [2] also shows that most of the heat conductivity measure- 
ments fall within a band of hx - ±0.10. 

One can recognize that the radiative heat-conductivity term 
is not negligible for temperatures above ~ 200 K. However, the 
first term on the rhs of Eq. |2] which describes the heat conduc- 
tivity through the particle network, still contributes consider- 
ably to the heat conductivity even for the highest temperatures 
in Fig. |2]so that a calibration of the heat conductivity model is 
possible. 

We also allowed a variation of the scaling parameter ei in 
Eq. |3] and tested whether the agreement between the thermal 
conductivity measurements of lunar regolith and the model can 
be further improved. We get a slightly better reduced chi- 
square if we vlss x - 0.37 ± 0.02 and ei - 1.66 + 0.21. As 
this agrees with the model with only one free parameter within 
~ 1.5 standard deviations, we henceforth used;^^ = 0.41 + 0.02 
andei = 1.34 + 0.01. 

5. Important properties of the analyzed objects 

Tab. [3] summarises important properties of the analyzed ob- 
jects, including the spectral classes of the objects, the diameters 
D of the objects, the gravitational acceleration g at the surface 
of the objects, the measured thermal inertiae F of their surface 
materials, their surface temperatures T during observations, and 
the estimated mean grain radii r of the surface regolith. The as- 
teroids Lutetia and Cybele can be assigned either to the spectral 



class C or M. In this work, these asteroids are treated as C-class 
objects. The diameters of the analyzed objects were either taken 
from literature or from the JPL Small-Body Databas^ 

The gravitational acceleration at the surface was calculated 
using the mass and the effective diameter of the object. If no 
measurement of the mass was available, we calculated the grav- 
itational acceleration assuming spherical objects with a mean 
density depending on the class of the asteroid: p — 1400 kg m~ 



(C), p = 2690 kg m-^ (S) and p = 4700 kg m"^ (M, Britt et al. 



2002 1. For Mercury the gravitational acceleration was taken 
from the NASA webpag^ 

For the Moon, the thermal inertia measurements were per- 
formed for temperatures between 200 K and 390 K. Thus, the 
uncertainty of the grain size estimation is dominated by the tem- 
perature variation. For Mercury, the influence of the tempera- 
ture (ranging from 350 K to 700 K) on the error of the particle 
size estimation was also taken into account. 

For objects with no direct measurement of the surface tem- 
perature, the maximum surface temperature was calculated us- 
ing the standard thermal model (see |Leboftky£tar 1986 and 
references therein for details). This model assumes spherical, 
slow rotating objects heated by sunlight and cooled by radia- 
tion. Each point on the surface is in instantaneous equilibrium 
with the Solar radiation and no heat is conducted through the re- 
golith towards the center of the object. This approximation was 
used due to the fact that all analyzed objects have long rotation 
periods (rotation period > 4 h). Following the standard thermal 
model, the maximum surface temperature at the subsolar point 
is given by 



T = 



cr erj 



(1 - A) 



1/4 



(8) 



( [Lebofsky and Spencer| [l989 |. Here, 5 = 1367 Wm^^ is the 
Solar constant, A is the Bond albedo of the object, is the 
semi-major axis of the Earth's orbit, a is the distance of the ob- 
ject to the Sun during observation, cr is the Stefan-Boltzmann 
constant, e is the emissivity of the surface material and rj is the 
so-called beaming parameter of the object, which takes devia- 
tions from the standard thermal model (e.g. non isotropic heat 
radiation) into account. In our work, we use t] - \ (spherical 
objects, standard thermal model). The bond albedo Ai,o„d was 
calculated using the phase integral q and the geometric albedo 
Ageom, Ahond - qAgeom- The values for the geometric albedo 
and the distance of the object to the Sun at the date of obser- 
vation were taken from the JPL Small-Body Database'. An 
approximation for the phase integral for asteroids is given by 
q - 0.290 + 0.684 G, with a slope parameter of G = 0.25 (see 
appendix in Bowell et al. 1989| l. Tab. [4] summarizes the used 
values and the derived temperatures. 

6. Determination of the regolith particle size 

Applying the thus validated heat-conductivity model to other 
small Solar System bodies, the mean size of the planetary re- 



'http://ssd.jpl.nasa.gov/sbdb.cgi 

-http://nssdc.gsfc.nasa.gov/planetary/factsheet/mercuryfact.html 
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Table 3: Important properties and estimated grain sizes of the analyzed objects. 
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Table[3](cont.): Important properties and estimated grain sizes of ttie analyzed objects 
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Table 4: Parameters used for computation of the surface temperature. 



Object 


4 




Date of Observation 


a [AU] 


r [K] 


1998 WT24 


0.56 


0.26 


2001-12-04 


1.01 


339+ 


1999 JU3 


0.06* 


0.03 


2007-05-16 


1.41 


277 


1996 FG3 


0.04* 


0.02 


2011-01-19 


1.38 


281 


Betulia 


0.08 


0.04 


2002-06-02 


1.14 


307 


Eros 


0.25 


0.12 


1998-06-29 


1.62 


252 


Elvira 


0.28 


0.13 


1983-07-28 


2.64 


197 


Bohlinia 


0.20 


0.09 


1983-08-09 


2.88 


191 


Unitas 


0.21 


0.10 


1983-07-31 


2.17 


219 


Dodona 


0.16 


0.07 


1983-07-11 


2.55 


226+ 


Nysa 


0.55 


0.25 


1983-07-27 


2.47 


218+ 


Thyra 


0.27 


0.12 


1983-04-28 


2.45 


205 


Pomona 


0.26 


0.12 


1983-07-31 


2.81 


191 


Lydia 


0.18 


0.08 


1983-06-25 


2.87 


213+ 


Ekard 


0.05 


0.02 


1983-06-13 


1.92 


238 


Herculina 


0.17 


0.08 


1975-01-01 


3.13 


183 


Cybele 


0.05 


0.02 


1985-12-10 


3.80 


169 


Vesta 


0.42 


0.19 


1970-01-01 


2.46 


200 


Pallas 


0.16 


0.07 


1970-01-01 


3.42 


176 



^^eom - geometric albedo. 
^hond- Bond albedo. 

a: distance to the Sun during observations. 
T: surface temperature at the subsolar point. 
^ : temperature calculatedj. 
* : geometric albedo fromU 



*: geometric albedo fromlLeon et al. 2011 




- 0.66 (M class asteroids). 
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golith particles can be estimated. This can be done by com- 
paring the heat conductivities derived from the thermal inertia 
measurements with those predicted by the model, treating the 
unknown volume filling factor of the material as a free parame- 
ter (see Sect. |2]i. 

To test this method, the grain size of the lunar regolith 
was estimated from thermal inertia measurements of the lu- 
nar surface (Fig. [3^). The dashed horizontal lines denote the 
heat conductivities derived from the measure d thermal iner- 
tia, r = 43 J K ' m-2s" 5 ( |Wesselink| |T948t , for filling fac- 
tors between (f> = 0.1 (uppermost line) and = 0.6 (low- 
est line). For the Moon, we derived the heat capacity c of 
the regolith material from the equation for the thermal iner- 
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Figure 3: Grain size estimation for tlie surface regolitli of tlie Moon and aster- 
oid (25143) Itokawa, Results of the regolith heat conductivity model (dotted 
curves) are shown for the lunar surface (a) and the surface of asteroid (25143) 
Itokawa (b) for different packing fractions of the regolith particles, = 0.1 
to (p = 0.6, in intervals of A0 = 0.1. The heat conductivities derive d from 
the thermal inertia measurements, r = 43 J K"' 
|T948landr = 750!^|]n J K"' m^^ s-^-^ (Itokawa; 



' m - s £f_(Moon; Wesselink| 



shown by the dashed lines, also for i 



0.1 tO( 



MilUer et al.| |2005 
■TO 



are 

Ihe restllting grain 



sizes are r = 48^2^ /.im for the Moon and r = 21;';|^ mm for asteroid (25143) 
Itokawa (crosses). Here, the size errors for the lunar regolith derive from the 
temperature variations of the lunar surface during the measurements, whereas 
the errors for (25 143) Itokawa stem from the uncertainties of the measured ther- 
mal inertia. For comparison, the 10% and 90% values of the size distribution of 
the lunar regolith samples returned with the Apollo- 1 1 mission ( McKay et al.| 
|2009 1 are also shown (dashed-dotted vertical lines). The model parameters used 
for the derivation of the regolith grain sizes are summarized in Tabs.^and|3] 



tia (see Eq. [TJ. Using the measured thermal inertia o f the 
lunai- surface, F = 43Jm 's JWessehnkl 1 19481, the 



measured heat conductivity of the lunar samples at 300 K, 
A - 2 X 10"^ Wm"' K"', and the values for the material and 
packing density of the lunar regolith shown in Tab. |2j we ar- 
rive at c = 711 J kg ' K ' . For the other analyzed objects, the 
values of the heat capacity summarized in Tab. [T] were used. 
All relevant values for the estimation of the grain sizes are pre- 
sented in Tabs. [T]and|3] 

The dotted curves in Fig. |3^ represent the model (Eqs. |2]|5]l 
for different volume filling factors, <p - Q.l to (p = 0.6. The 
intersection between the model curve and the respective hori- 
zontal line yields the mean regolith particle size for each as- 
sumed packing density (crosses in Fig. [3^). As can be seen 
from Fig. |3^, one can find such an agreement only for the lower 
three packing fractions. Taking into account the uncertainties 
in measuring the thermal inertia of the lunar regolith, we thus 
arrive at an average grain size of r = 48^2^^ /vm. This falls well 
into the measured size distribution of the lunar regolith samples 
returned with the Apollo-11 mission (McKay et al. 2009 [ l. In 
Fig. [3^ we show the 10% and 90% values of the size distribu- 
tion of the lunar regolith by the dashed-dotted vertical lines. 

As can be seen, there is no agreement between the model and 
the derived heat conductivities from thermal inertia measure- 
ments for packing fractions above (p = 0.3. This is in contra- 
diction to the measurements of the volume filling factor of the 

[T972l l. 



lunar regolith, which indicate (p ^ 0.5 (Mitchell et al. 
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Figure 4: Grain size estimation peii'ormed for the surface regolith of the Moon 
using a constant value for the volume filling factor of (p = 0.5. We varied the 
thermal inertia in order to find a match between the corresponding heat conduc- 
tivity and the heat conductivity resulting from the model. For a thermal inertia 
of r = 60 J m"^ K"' s""-^, a match between the measured heat conductivity, 
A = 2.7 X 10"^ Wm"' K"' (dashed line) and the modeled heat conductivity 
(dotted curve) was found. The comparison of the heat conductivities yields a 
mean grain radius of r = 44fim. For comparison, the 10% and 90% values 
of the size distribution of the lunar regolith samples returned with the ApoUo- 
11 mission (McKay et al. 2009 1 are also shown (dashed-dotted vertical lines). 
Additionally, the heat conductivity estimated from the measurement of the in- 
frared radiation of the lunar surface during the Apollo- 17 mission is shown, 
A = 1.5 X 10^3 Wm^'K"' (da.shed-dotted-dotted line). 
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Thus, we increased the thermal inertia in order to find a match 
between the derived heat conductivity and the heat conductivity 
resuhing from the model (see Fig. |4]). We only get an agreement 
if we use a thermal inertia of F = 60Jm"^K"' s"*^^, which 
results in a heat conductivity of A = 2.7 x lO^-'Wm^'K"' 
(dashed line in Fig. [4]i. We thus obtained a regolith grain ra- 
dius of r = 44 ij.m. The derived grain size is very close to the 
grain size estimated by varying the volume filling factor of the 
material, but having volume filling factors lower than cp - 0.5. 
However, our value for the heat conductivity is almost a factor 
of two higher than the value derived from measuring the in- 
frared radiation of the lunar surface during the Apollo- 17 mis- 
sion, A - 1.5 X 10 Wm ' K ' (dashed-dotted-dotted fine in 
Fig. |4) [Keihm et al] [T973| ). [Keihm et al^ ( [T973] l derived the 
heat conductivity from a nighttime cool-down curve starting at 
~ 130K and decreasing to ~ 105 K. Thus, the discrepancy 
between our result (derived for T - 295 K) and the value esti- 



mated by Keihm et al. ( 1973 1 is caused by the usage of different 
temperatures as can be seen in Fig. |2] 

After these calibration and validation steps, we applied our 
model to various objects for which thermal inertia measure- 
ments were available (see Tab. |3]l and estimated the grain size of 
the surface regolith. As the volume filling factor of the regolith 
is a priori unknown, we treated its value as a free parameter 
in this work. Fig. [SJ) exemplarily shows our analysis for the 
surface regohth of the asteroid (25143) Itokawa, target of the 
sample return mission of the Hayabusa spacecraft. For the mea- 

( ,Muller, 



sured thermal inertia of F = 750;^3qq J K ' m ^ s "-^ 



et al. 



2005|l, we arrive at a mean grain radius of r = 21" 



which is in agreement with in-situ observations (Yano et al. 



14 mm. 



2006 Kitazato et al. 2008 1. The uncertainty of the grain size 



measurement is completely dominated by the error of the ther- 
mal inertia determination. In contrast to the Moon, the heat 
transport in the regolith on asteroid (25143) Itokawa is gov- 
erned by radiation so that uncertainties in the material proper- 
ties are unimportant. 

Figs. |5]|8]show the grain size determination for all other ob- 
jects listed in Tab. |3] The surface temperature of the objects 
was taken from the literature or was calculated (see Sect. [5]l if 
no direct temperature measurement was available. 

For the Moon and Mercury, the thermal inertia measurements 
were performed for a wide range of temperatures (see Tab. [3]l. 
For the Moon, we calculated the mean particle size using a sur- 
face temperature of T = 295 K. For Mercury, a surface tem- 
perature of r = 700 K was used for the estimation of the mean 
grain size. In these two cases, the temperature variation causes 
an additional uncertainty of the grain size estimation. Both er- 
rors stemming from the temperature range are on the same order 
of magnitude. Hence, two error bars are later shown in Fig. [9] 
for Mercury, one for the uncertainty of the grain size estimation 
due to the temperature variation and the other for the uncer- 
tainty due the error of the thermal inertia measurement. For the 
Moon, unfortunately no error discussion of the thermal inertia 
estimation is available (Wesselink 1948). Thus, the uncertainty 
of the lunar regolith grain size estimation is determined by the 
temperature range only. 

In some cases, the grain size estimation was not possi- 



ble for the entire error interval of the measured thermal iner- 
tiae. For example, the thermal inertia of th e asteroid Dodona 
was estimated to be 83 



68 J m-^ K 



(Delbo and Tanga 



|2009) , and we only found solutions for thermal inertiae ranging 
from36Jm"^ K"^ s""^ to 151Jm"^ K"' s"°^. For values below 
36 J m"^ K"' s"°^, the required values for the density of the sur- 
face regolith are too high. Tab. |5] summarizes the cases, where 
only part of the measured thermal inertia range could be used 
to estimate the grain size of the regolith. The range of thermal 
inertiae used in our calculations is represented by F' in Tab. |5] 



7. Functional dependence of the mean regolith grain size on 
the gravitational acceleration of the asteroid 

We applied our model to a number of asteroids for which 
thermal inertia measurements were available (see Sect. |6]l. 



Delbo and colleagues (Delbo et al. 2007 Delbo and Tanga 



|2009) compiled such results and found that an inverse correla- 
tion between the thermal inertia and the diameter of the asteroid 
exists. We used these data and augmented them by thermal in- 
ertia measurements of the asteroids 1998 WT24, 1999 JU3 (the 
possible target of the Hayabusa 2 mission), 1996 FG3 (the pos- 
sible target of the Marco Polo-R mission) and (2867) Steins as 
well as of the Martian moons Phobos and Deimos and the planet 
Mercury. We sorted these objects into classes, with each class 
having its own set of material parameters, i.e. a temperature 
dependent heat conductivity and material density (see Tabs. [T] 
and[3]l. As done for the Moon, we varied the unknown packing 
density of the regolith from (f> - Q.I to (p - 0.6 and thus derived 
six diff'erent heat conductivity values for each body. We then 
applied our model (Eqs. [2]|4]i and searched for intersections be- 
tween the model and the respective heat conductivities, which 
yielded the grain size of the regolith (see Figs. [3]and[5][8]). We 
calculated the mean grain size by arithmetically averaging the 
(up to six) solutions. In the rare cases of double intersections, 
we chose the larger of the two possible values of the grain size. 
In Fig. |9] we show the results of our analysis. As for the thermal 
inertiae, the derived mean particle size of the regolith follows 
an anti-correlation with the diameter and, thus, with the grav- 
itational acceleration of the body (see Tab. [3]). Small objects, 
like Asteroid (25 143) Itokawa with a diameter of only 320 m 
and a gravitational acceleration of only ~ 9.3 x 10"^ m s"^, pos- 
sess a much coarser regolith, with typical particle sizes in the 
millimeter to centimeter regime, whereas those bodies exceed- 
ing 100 km in size {g ^ 5 x 10"^ ms"^) possess a much finer 
regolith, with grain sizes between 10 fim and 100 yum. The un- 
certainties in the determination of the particle size (error bars in 
Fig. |9| are dominated by the errors of the thermal inertia mea- 
surements, which are typically 50% of the mean value (see Tab. 

The decrease in grain size can qualitatively be understood as 
a result of the collision history of the asteroids. As mentioned 
above, previous studies have shown that hyper-velocity impacts 
lead to a size discrimination of the regolith, which consists of 
those impact fragments whose velocities did not exceed the es- 



cape speed of the parent body (Chapman 1976 Housen and 
Wilkening 1982 Chapman 2004[ l. As naturally the smaller 
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Figure 5: Grain size estimation for tlie surface regolitli of 1998 WT24 (a), 1999 JU3 (b), 1996 FG3 (c), Steins (d), Betulia (e) and Deimos (f). 
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Figure 6; Grain size estimation for the surface regolith of Eros (a), Phobos (b), Elvira (c), Bohlinia (d), Unitas (e) and Dodona (f). 
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Figure 7: Grain size estimation for the surface regolitli of Nysa (a), Tliyra (b), Pomona (c), Lydia (d), Ekard (e) and Lutetia (f). 
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Figure 8; Grain size estimation for tlie surface regolitli of Herculina (a), Cybele (b), Vesta (c), Pallas (d), Ceres (e) and Mercury (f). 
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Table 5: Thermal-inertia ranges measured (F) and used in this work (T') 



Object 

Dodona 

Herculina 

Cybele 

Pallas 



-2 ^-1 ,-0.5n 



r[Jm-^K-' s 

-68 

[5+8 
'-^ -8 

15- 
10- 



jDelbo and Tanga 2009 



Miiller and Lagerros 


1998 


Miiller and Blommaert 


2004 


Miiller and LagerrosI 


1998 



r [Jm 

38 - 151 
12 - 23 
9-23 
10 - 15 



F: measured thermal inertia, 
r' : used thermal inertia. 



fragments gain the highest velocities in impacts (Fujiwara and| 



Tsukamoto 1980[ [N akamura and Fujiwara 199 



Nakamura 



1993 1 [Nakamura et al.| |1994) [Vickeryj [19861 [1987) and as the 
asteroids were roughly subject to the same impactor sizes and 
velocities, the expectation that small bodies are covered by the 
coarser fragments, whereas large bodies were able to hold on 
to most of the size distribution, seems to be justified (it should, 
however, be borne in mind that electrostatic effects might also 
play a role for the release or capture of surface dust (jLeeJ 
1996[ l). Typical impact generated size distribution functions 
have the tendency to be dominated in number by the smallest 
fragments ( [Nakamura et al.| |1994^ so that the mean grain size 
responsible for the thermal conductivity of the regolith is ex- 
pected to be close to the smaller limit. 

To analyze the data, we fitted the logarithmic data with the 
following functions. 



log 
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and 



log 



1 yum 



1 + exp[/72(log[j^] - ^3)] 



+ bA (11) 



(Figs. [9^-c). Tab. [6| summarizes the used fit functions, the val- 
ues of the fit parameters, and the resulting values of the reduced 
chi square, x^^.^- The power law fits ri oc g^^l^ and r2 oc g^'^^- 
are motivated by the fact that laboratory experiments (Fujiwara 
and Tsukamoto' '1980'; 'Nakamura an d Fujiwara[ [1991( |Naka-| 
mura, 1993, Nakamura e t al. 1994) and studies of the crater 
structures on the Moon ( |Vickery[ |1986 19^7]l have shown that 
hyper- velocity impacts accelerate small fragments to higher ve- 
locities than large ones. Fig. 5 in Nakamura et al. ( 1994[ l 
shows a comparison of the laboratory data and the results of 
the crater structure analysis. A simple interpolation of the data 



on the upper left and lower right of Fig. 5 in Nakamura et al. 



( 1994 1 by a power law function yields a dependency between 
the velocity Vfrag and the size dfr ag of the ejected f ragments of 



df^'^. Note that Fig. 5 in 



Vfrag - "frag ■ 

alizes the velocity of the fragments on tne vertical axis and the 
size of the fragments on the horizontal axis. However, our Fig. 



Nakamura et al. 



(19941 visu 



[9]shows the size of the regolith on the vertical axis and the grav- 
itational acceleration, which is proportional to the escape speed 
of the body, on the horizontal axis. We have inverted the esti- 
mated exponent in order to apply the power law function to our 
data. Hence, the grain size of the asteroidal regolith particles 
should follow a power law function of the surface acceleration 

with an exponent of 2/3. The smoothed step function r^ig) 

was chosen because of the appearance of the data, which sug- 
gest a rather steep decrease of the regolith grain size for bod- 
ies of ~ 100 km in size {g -~ 5 x 10"^ ms"^), close to the 



value above which asteroids are thought to be primordial (Bot- 
|tkeetal.||2005] l. 

All fit functions were applied to two different data sets (DS 1 
and DS2). The first data set (DSl; filled circles in Fig. [9| 
consists of all asteroids without the asteroids analyzed using 
the IRAS data (short: IRAS asteroids). The second data set 
(DS2; filled circles plus crosses in Fig. [9|l consist of all aster- 
oids including the IRAS asteroids. The binary asteroid 1996 
FG3 (bowtie), the Martian moons Phobos and Deimos (filled 
squares), the Moon (downward triangle) and Mercury (upward 
triangle) were excluded from the fits because the surface re- 
golith on these objects was probably influenced by the presence 
of the close-by asteroid (1996 FG3) or planet (Mars or Earth), 
or by the location deep within the solar gravitational potential 
(Mercury). 

We found that fixing the exponent of the power law to -2/3 
results in a slightly better fit than using the exponent as a free 
parameter (see Tab. [6]|. On top of that, both power laws are 
identical within the errors of the fit parameters. As was already 
addressed in [Delbo and Tang a (2009), the thermal inertiae of 
the asteroids observed by IRAS seem to be offset and to pos- 
sess a steeper slope than the rest of the sample. Excluding this 
sample from the above analysis yields a better fit of the power 
law functions to the data. 

However, the smoothed step function yields a better fit to 
the data if the asteroids observed by IRAS are taken into ac- 
count (see Tab. [6]). Excluding the IRAS asteroids decreases the 
quality of the fit in case of the smoothed step function. The 
residuals of the smoothed step function including the IRAS as- 
teroids are almost identical to the residuals of the two power 
law functions without the IRAS data. The smoothed step func- 
tion was chosen due to the appearance of the data and has no 
physical explanation yet. However, electrostatic effects close to 
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Figure 9: Dependence of the regolith grain size on tlie gravitational acceleration of tlie asteroid. The derived grain sizes for the regolith of the asteroids in the 
sample of|pelbo et al. (2007J and |Delbo and Ta nga (20091 as well as for the asteroids 1998 WT24, 1999 JU3 (the possible target of the Hyabusa 2 mission), 1996 
FG3 (the possible target of the Marco Polo-R mission), (2867) Steins, the Martian moons Phobos and Deimos, and the planet Mercury are shown as a function 
of the gravitational acceleration of the planetary body. Three different fit functions (dotted and dashed curves), ri(g) (a), r2ig) (b) and rjig) (c), were fit to two 
different data sets (see Tab. [6|. The first data set (DSl) consists of all asteroids without the IRAS asteroids (filled circles). The second data set (DS2) consists of all 
asteroids including the IRAS asteroids (filled circles plus crosses). The binary asteroid 1996 FG3 (bowtie), the Martian moons Phobos and Deimos (filled squares), 
the Moon (downward triangle), and Mercury (upward triangle) were excluded from the fits. The green colored data denotes the measured size distribution of the 
lunar regolith samples ( McKay et al. 2009 1 and the arrow denotes the lower limit of the grain radius on the surface of Itokawa jYano et aL]|2006[|Kitazato et al.| 
[2008^ . Additionally, the residuals of the data to the fit functions are shown in logarithmic units for DS 1 (lower row) and DS2 (upper row). 



Table 6: Fit functions and resulting fit parameters used to describe the anti-correlation between regolith particle size and gravitational acceleration of the asteroid. 



Fit Function 


Data Set 


b, 


b2 


63 






Eq. 
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DSl 


1.23 ± 0.14 








0.24 


Eq. 
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DS2 


1.61 ± 0.14 








0.39 


Eq. 
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DSl 


1.08 ± 0.27 


0.74 ±0.12 






0.26 


Eq. 


10 




DS2 


1.60 ± 0.30 


0.67 ± 0.14 






0.42 


Eq. 


11 




DSl 


4.43 ± 0.00+ 


0.75 ± 0.00* 


-2.17 ± 0.00* 


0.48 ± 0.00* 


0.32 


Eq. 


11 




DS2 


2.30 ± 0.00* 


2.90 ± 0.00* 


-1.26 ± 0.00* 


1.31 ± 0.00* 


0.31 



^red- reduced chi square. 

DSl: data set 1; all asteroids excluding the IRAS asteroids. 
DS2: data set 2; all asteroids including the IRAS asteroids. 

^: the unceitainties of the fit parameters are: Ab^ = l.OOx lO-"*, Ab2 = 4.94 X 10"^, A/?3 = 7.63 x 10"'^ andA;74 - 5.20 x IQ-'^. 
^: the uncertainties of the fit parameters are: Abi = 5.27 X 10"^, Abj = 3.15 X 10""^, A63 = 3.56 x 10"^ and Ab4 = 3.82 x 10"^ 
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the surface of the small asteroids have been mentioned to influ- 



ence the escape of ejecta ( Lee 1996 1. The smoothed step func 



tion obviously shows a separation of the objects into two dif- 
ferent classes. Small objects with diameters less than ~ 20km 
are covered by relatively coarse regolith 
grains with typical particle sizes in the millimeter to centime- 
ter regime, whereas large objects with diameters in excess of 
~ 80 km ^ 6 X 10"^ ms"^) are covered by very small re- 
golith particles with grain sizes between 10 fim and 100 fim. 



8. Influence of the network heat conductivity on the grain 
size estimation 



We determined the grain size of planetary regolith for a large 
number of asteroids, the Moon, the Martian moons and Mer- 
cury with diameters between 0.3 km and 4,880 km (see Tab. [3]l 
and find an anti-correlation between the regolith grain size and 
the gravitational acceleration of the planetary body (see Fig. |9]l. 
The anti-correlation between grain size and gravitational accel- 
eration is independent of the spectral class of the asteroids (see 
Fig. [TT] i. This outcome supports the idea that planetary regolith 
is formed by hyper-velocity impacts, which have ground down 
the initially rocky material to ever finer particle si zes (|Chap- 



man| [T976i [Housen and Wilkening| [T982t |Chapman| [2004 1. 
However, the determined grain sizes also show a separation 



As can be seen above, small asteroids (< 20 km) are cov- 
ered by relatively coarse regolith (in the millimeter to cen- 
timeter regime, see Sect. |7]i. In the case of coarse regolith, 
the total heat conductivity of the regolith is dominated by ra- 
diative heat transport. Here, we show that the regolith grain 
size of small asteroids (< 50 km) can be easily estimated by 
only calculating the radiative heat transport of the regolith, 
A X /liad = 8crer^ A(r, 0) (see Eq. This approximation 
provides grain sizes with an accuracy of approximately a factor 
of 2. 

Fig. 



10 shows the normalized grain radius r' /r as a func- 



tion of the diameter of the object. Here, r' is the derived grain 
size of the surface regolith particles under the assumption that 
the heat transport is due to radiation, i.e. /lsoiid(7^) H(r, T, (p) <& 
^creT^ A(r, (p). If the heat conductivity was solely determined 
by radiative transport, all data would be located on the dashed 
line, r' /r - 1. The deviation of the data from the dashed line 
shows the influence of the heat conductivity of the solid net- 
work of regolith on the grain size estimation. For small ob- 
jects (< 50 km), the deviation of the relative grain radius to the 
dashed line is relatively small (less than a factor of 2). This 
demonstrates that for small objects (< 50 km), the total heat 
conductivity of the regolith is generally dominated by radia- 
tive heat transport. In this case the estimated grain sizes are 
independent of the specific material properties, as long as the 
condition e s; 1 still holds (see Eqs. |2]and[5]l. For larger ob- 
jects (> 50 km), the heat conductivity of the solid network of 
regolith and, thus, the material properties are not negligible. 



9. Conclusion and discussion 

In this work, we presented a new method to determine the 
grain size of regolith particles by using remote measurements 
only. We utilized measurements of the thermal inertia per- 
formed for the regolith of various objects in the Solar System 
(see Tab. [3]). With the knowledge of the thermal inertia, the 
surface-material properties and assumptions about the packing 
density of the regoUth particles, the heat conductivity of the 
surface regolith was derived. The grain size of the planetary 
regolith was then determined from a comparison of the derived 
heat conductivity with a modeled heat conductivity of granular 
materials in vacuum. 
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Figure 10: This figure shows the influence of the solid network heat conductiv- 
ity on the grain size estimation by comparing the normalized grain radius r' /r 
with the diameter of the object. The grain radius r' was estimated under the 
assumption that heat is only transported by radiation, i.e. /isoiidC^) H(r, T, <j>) <s. 
ScreT^ A(r, <p). The deviation of the data from the dashed line shows the im- 
portance of the heat conductivity of the solid network of regolith on the grain 
size estimation. 
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Figure 1 1 : Comparison of the derived grain sizes for the different spec- 
tral classes of the asteroids: S-class (pluses), C-class (squares) and M-class 
(crosses). 
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of the objects into two different regimes. Small objects with 
diameters less than ~ 20 km < 1 x 10"^ ms"^) are covered 
by relatively coarse regolith grains with typical particle sizes 
in the millimeter to centimeter regime and large objects with 
diameters bigger than ~ 80 km ig ^ 6 x 10"^ ms"^) possess 
very small regolith particles with grain sizes between 10 fim and 
100 fim. A previous estimation of the grain size of asteroidal 
surface regolith derived from polarization measurements (grain 
size: ~ 100/zm, Dollfus andZellner 1979| l is in agreement with 
our result. 

Our method to determine the grain size of planetary regolith 
depends on the measured thermal inertia of the object, the sur- 
face temperature of the body during observation and the mate- 
rial properties of the regolith. Due to Eq. [T] the uncertainties of 
the thermal inertia measurements, which are typically 50 % of 
the measured value, primarily determine the error of the grain 
size estimation (see Tab. [3]l. Typical temperature variations of 
the analyzed objects have only a minor influence on the grain 
size estimation (less than ~ 10% of the estimated value). Only 
for the Moon and Mercury, the error of the grain size estimation 
due to the temperature variation becomes important (see Sect. 
|6]l. The used material properties only have an influence on the 
grain size estimation if the diameter of the object is larger than 



50 km (see Fig. 10 in Sect, p 



In principle, a rough estimation (accuracy within a factor of 
two) of the grain size of the surface regolith of small objects 
can be performed by assuming radiative transport only (see Fig. 
10 1. For bigger objects and, thus, smaller regolith particles, the 
heat conductivity of the regolith network becomes more impor- 
tant so that the uncertainty of the grain-size estimation stem- 
ming from this simplification increases. 

We conclude that our method allows the direct determination 
of the grain size of asteroidal or planetary regolith from mea- 
surements of the thermal inertia of the parent body. The derived 
mean particle radii vary from millimeters to centimeters for km- 
sized objects to 10 - 100 jum for bodies exceeding ~ 80 km in 
diameter This finding is important for the planning of future 
landing and/or sample-return missions as it allows narrowing 
down the physical state of the regolith by state-of-the-art re- 
mote observations. Smaller particle sizes in the regolith lead to 
a higher cohesion of the uppermost layers, which is, for small 
bodies, not determined by the local gravity field but dominated 
by inter-particle forces. 

On top of that, our method can be used to validate models 
of the formation and evolution of planetary regolith. Hyper- 
velocity impacts cause the ejection of impact fragments with 
a characterisitc size distribution and size-velocity relation ( |Fu-| 
jiwara and Tsukamoto' '1980; 'Nakamura and Fujiwara' '1991* 
Nakamura , 1993 ; Nakamura et al. , 1994;yickery, 1986 , 1987) . 



In combination with the impact geometry and the escape speed 
of the target object, one can thus derive a zero-age regolith size 
distribution. Space wheathering, granular transport eff'ects, and 
other internal and external influences might alter this size dis- 
tribution over time. 
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